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Abstract 

I investigate a new idea of perturbation theory in covariant canonical 
quantization. I present preliminary results for a toy model of a har- 
monic oscillator with a quartic perturbation, and show that this method 
reproduces the quantized spectrum of standard quantum theory. This re- 
sult indicates that when the exact solutions to classical equations are not 
known, covariant canonical quantization via perturbation theory could 
be a viable approximation scheme for finding observables, and suggests 
a physically interesting way of extending the scope of covariant canoni- 
cal quantization in quantum gravity 
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1 introduction 



There is a long-standing debate in quantum gravity on the merits of covariant versus canon- 
ical approaches to quantization. Canonical approaches require a split of spacetime into 
space and time, contrary to the spirit of general covariance. covariant approaches based 
for instance on path integral formulation, is intrinsically perturbative and require a clas- 
sical background. There is one approach, covariant canonical quantization, that avoids 
both these weaknesses, and provides a resolution to the covariant versus canonical debate, 
at least in principle. In this approach one quantizes the space of solutions to the classical 
equations. This solution space is isomorphic to the space of initial data (provided the initial 
value problem is well posed), and each initial datum is in turn a point on the phase space. 
Covariant phase space [1,2] provides a starting point based on observables of the theory, 
since canonical data are themselves observables in a time reparametrization invariant the- 
ory. In spite of this nice feature the method remains relatively unexplored, mainly because 
it relies on the knowledge of the full set of solutions to the equations of motion. To address 
this difficulty, it has been suggested [3] that it may be useful to extend the scope of this 
idea through perturbation theory. In this paper, I attempt such a construction. In particular 
I see if one can quantize a harmonic oscillator with a quartic perturbation by quantizing 
the space of solutions generated at each order via a classical perturbation expansion. 

General relativity differs from other classical field theories in many aspects. The clas- 
sical theory is nonlinear and its underlying symmetry, diffeomorphism invariance, implies 
that observables of the theory are non-local, a local definition of a field defined at a space- 
time point x M is not physically meaningful, since diffemorphism actively moves one space- 
time point to another. Diffeomorphism invariance also implies that there is no preferred 
background time. This raises problems for canonical quantization, which in its usual for- 
mulation requires a preferred splitting of spacetime into space and time. The absence of a 
preferred time makes the description of observables and dynamics one of the central chal- 
lenges in quantum gravity [3-9]. In the canonical formulation, one works on a fixed time 
slice. Here we usually start with the canonical geometrodynamic variables (q a b, n ab ) [10] 
or the Ashtekar variables [11]. The challenge then is to describe the dynamics of spacetime 
in an invariant way, after having chosen to work in the first place with an "artificial" split of 
spacetime into space and time. In classical general relativity the evolution from an initial 
time slice to a final one is independent of the choice of time in the intervening spacetime. 
But it is not clear if the same is true in the quantum theory — Torre and Varadarajan [12] 
recently showed that even in the absence of gravity, quantum evolution of a scalar field can 
be slice-dependent. 

In the covariant picture, we work within a path-integral framework and employ a back- 
ground field method that closely follows the standard approach to quantum field theories on 
a fixed background. Although the path-integral formally yields quantities that are covariant 
under diffeomorphisms of the background metric, the metric variables that we started are 
by themselves not diffeomorphism invariant. In contrast to the usual field theory approach, 
where correlators of locally defined fields make sense, in this case the correlators are com- 
plicated non-local functions of curvatures. If the perturbation h^ v is large enough, the 



1 



signature of the metric could change, and that will mean that the Lorentzian path-integral 
could miss the space of physical metrics. Further, Goroff and Sagnotti [13] have shown 
that gravity is nonrenormalizable at two loops — a result that effectively ended hopes of 
a conventional quantum field theoretical approach to quantum gravity. Conventional per- 
turbative evaluation of the path-integral thus seems untenable for quantum gravity. Hence, 
there is motivation to formulate a new approach to perturbation theory for quantum gravity. 

In covariant canonical quantization, also known as covariant phase space quantization, 
the starting point is the space of classical solutions. For a theory with a well-posed ini- 
tial value problem, the classical solutions are uniquely determined once arbitrary initial 
data is specified. Each set of initial data is also a point on phase space, and hence by 
restricting the classical solutions to an arbitrary initial value surface E , one achieves an 
isomorphism between the space of solutions § and the phase space T*Q (see section 2 
for details). The challenging aspect here is to find interesting physical observables — since 
the canonical data (q,p) are time independent, the wave functions ip(q) or (j>(p) in covari- 
ant canonical quantization are Heisenberg picture wave functions, and nothing seems to 
evolve. However, Rovelli [14-16] has shown that the observables in this framework are 
"evolving constants" — one parameter families of phase space functions with vanishing 
total "time" derivatives. As we discuss in section 2, the canonical data are themselves 
"evolving constants" in the sense defined by Rovelli. Time evolution is now encoded as 
a mapping between data. As noted above, the main obstacle in making progress with this 
idea (except in a few cases such as Carlip (2+1) [17, 18] and Torre [19] (cylindrical waves)) 
is that in general relativity in 3+1 dimensions we do not know the exact general solutions, 
hence motivating a perturbative approach to extend the scope of covariant canonical quan- 
tization. 

The proposal of the present paper can be outlined in the following way. I employ a 
classical perturbation expansion around a known solution to construct the canonical data 
of a different problem that cannot be solved exactly. Formally, this amounts to solving the 
Hamilton-Jacobi equations at each order, and then, in the spirit of the covariant canonical 
approach, identify the solutions at each order with canonical data of the full problem at that 
order. We can then pass onto a quantum theory by representing the perturbed canonical data 
as self-adjoint operators on the Hilbert space of basis states of the exactly soluble problem. 

Constructing observables by such "systematic approximations" is not altogether a new 
idea. It was first proposed in the context of general relativity by Bergmann et al [20]. Since 
then, the idea have been reiterated as an interesting alternative approach to the problem of 
observables by Carlip [3], Torre [9] and in a different guise by Gambini and Pullin [21]. 

However, covariant canonical methods, to my knowledge, have never been tested in 
perturbation theory. Given the success of this approach for the classical theory on one 
hand, and skepticism regarding its implications in quantum theory on the other (see com- 
ments in Ashtekar et al. [1]), testing this idea on a toy model as a preliminary investigation 
merits attention. As a first step, we construct the quantum theory of the quartic oscillator 
via classical perturbation theory. Section 3 of this paper carries the details of this construc- 
tion. Further,as we shall discuss in the concluding section, the quartic oscillator model may 
provide signposts on how the operational procedure can be adapted to derive meaningful 
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results for minisuperspace Bianchi models. Because of the high degree of symmetry in 
these models, they are effectively quantum mechanical problems, and provide the simplest 
nontrivial examples of dynamical spacetimes. There is also a large body of literature de- 
voted to reducing the Einstein equations to equivalent Hamiltonian dynamical systems for 
these models. It may be possible to exploit these results if we know the steps in quantizing 
a Hamiltonian system via the proposed approach. 



2 Quantum mechanics with classical canonical data 

In the simplest route to canonical quantization one starts with an underlying classical struc- 
ture and proceeds to define the Heisenberg algebra EI of quantum observables via a tran- 
sition from the Poisson bracket or symplectic 2-form * f2 = dp A dq to the commutator 
[q,p\. As emphasized by Crnkovic et. al [2], this process is not tied down to the seem- 
ingly non-covariant a priori choice of coordinates and momenta. Rather, the fundamental 
geometric structure on phase space, the symplectic 2-form, takes precedence. This follows 
from the observation that the presymplectic 1-form = p a dq a is a completely co variant 
structure that can be read off from the boundary term in the variation of the action when 
the variation is restricted to the classical solutions to the equations of motion [22,23]. In 
other words, when restricted to the classical solutions, 

f d tf 

^boundary = / dt—(p a dq a ) = [9] 



dt 

VL: = dQ (2.1) 

When the initial value formulation is well posed, classical solutions (q(t),p(t)) are gener- 
ated uniquely from initial data (q,p) by a canonical transformation C 

C : (q,p) !-»■ q{t) := q(q,p;t); p(t) := p(q,p;t) (2.2) 

By restricting the solutions to an arbitrary initial value surface E, this canonical transfor- 
mation can then be inverted to represent data in terms of solutions, that is we can find 
expressions, q = q(q,p; t),p = p(q,p; t). Hence, the space of solutions of the equations of 
a motion is isomorphic to the space of initial data. However, each of the data points is also 
a point on the phase space. Thus the symplectic 2-form on phase space Cl can be pulled 
back to the space of data § to define a symplectic 2-form co = dq A dp = C*f2 on the latter. 
Once this is done, quantization can be performed by a finding suitable self-adjoint opera- 
tors q and p, which generates the Heisenberg group H. Wave functions ip(q) or cf>(p) are 
defined for all times and are thus the wave functions of the Heisenberg picture of quantum 
mechanics. In passing let us also observe that, since data is independent of time, the total 
time derivative vanishes, that is 



*See Isham [28] for a review of algebraic quantization with the symplectic structure as a fundamental 
building block. 
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An "evolving constant"^, introduced by Rovelli [14] is a one parameter family of func- 
tions, parametrized by "time" t that satisfies .4 = 0. Hence from (2.4) we see that canonical 
data satisfies this criterion and hence are themselves evolving constants. 

When the full classical solutions are not known, we can re-formulate the problem of 
finding solutions in a perturbative sense outlined below. To be concrete we will confine 
ourselves to a system with a Hamiltonian description. Let us assume that we have a system 
described by a Hamiltonian H for which the full classical solutions are known. t This in 
particular means that the observables of the unperturbed problem (qo,Po) are known and 
that the unperturbed symplectic structure ujq = dpo A dqo is given. Assume then a small 
perturbation Hi, so that the full Hamiltonian is 

H = H + eH x (2.4) 

The idea is to use the full Hamiltonian to "evolve" the dynamical system, and in particular 
construct the observables of the problem described by H in terms of the canonical data of 
the unperturbed problem. The canonical variables of the full problem are then obtained as 
a perturbation expansion defined as, 

q = q + eq 1 + e 2 q 2 + .... 

p = p + epi + e 2 p 2 + .... (2.5) 

From the Hamilton's equations and the requirement that the observables of the full theory 

satisfy < q,p > = 1 at each order, we obtain solutions for (q x , q 2 ) and (vi,p2) and con- 
l J PB 

struct the solutions defined above at each order as functions of (go, Po) . and investigate the 
consequences of the perturbation on the zeroth order symplectic structure. The latter when 
promoted to a commutator algebra for the corresponding operators, encode the "evolution" 
of the unperturbed basis states that define the Hilbert space on which the quantum theory 
is defined. We now turn to the details of application of this scheme to the quantum theory 
of the quartic oscillator. 



3 Quantum theory of quartic oscillator via classical per- 
turbation 

The complete Hamiltonian for the toy model is described by 

2 2 2 4 

H = P - + ^ + ^ (3.1) 
2 2 4 

The basic canonical variables of the unperturbed harmonic oscillator problem will be de- 
noted by (qo,po). The classical solutions of the unperturbed Harmonic oscillator is de- 
noted as q(t) := q cosut + ^ smut. In particular this means that for t = 0, q = q and 

^In the proposal put forward by Bergmann et. al [8], it is assumed that the lowest order solution is the 
Minkowski metric, and the expansion is "in terms of a deviation from this trivial constant solution". 
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p = Po represent oscillator data. Time dependent solutions will be denoted with an overbar 
through out the remainder of this text, while their restriction to the t — surface will be 
denoted without it. The steps leading to the quantum theory can be outlined as follows: 



1. Hamilton's equations, in conjunction with the perturbation expansion 



q{t) = q (t)+eq 1 (t) + e 2 q 2 (t) + ... 
p{t) = p {t) + epi(t) + e 2 p 2 {t) + ... 



(3.2) 



are used to construct the solutions as q = q(qo,Po; t),p = p(qo,Po', t). 

2. We impose {q(t) , p(t)}p B = 1 to all orders in perturbation theory, and use this con- 
dition to obtain an expression for the Poisson bracket of the oscillator variables. We 
find that the Poisson algebra is preserved under time evolution, hence we work with 
physical quantities defined at t = 0. The perturbation manifests itself as a deforma- 
tion of the classical Poisson algebra of (qo,po). The Hamiltonian (3.6) is expanded 
to each order in perturbation theory with the aid of the solutions. The equations of 
motion ensure that it is an evolving constant in the sense of section 2. 

3. We define a set of new canonical variables, expressed as functions of (qo,Po) at each 
order in perturbation, which satisfies a classical canonical algebra. The introduc- 
tion of this canonical pair is motivated by the adjustment of the deformation in the 
classical Poisson algebra of the oscillator variables (see section 3.1 for an explana- 
tion on the meaning of the term adjustment used here). This enables a transition to 
a canonical pair of operators that generates a Heisenberg group EI at each order of 
perturbation. This is needed to construct the Hilbert space of states and to facilitate 
operator ordering of the Hamiltonian. 

4. We promote the latter variables with the choice of Weyl ordering to hermitian op- 
erators acting on the Hilbert space H osc of the unperturbed oscillator states, and 
construct a Hamiltonian operator H whose action on these states gives the quantized 
spectra at each order of perturbation. 

Proceeding in step 1, Hamilton's equations lead to the differential equations 



These must be satisfied at each order in e, and in particular yield the equations of motion 

for for qupi, q~2, P2 respectively as 



p = —uj 2 q + eg 3 ; < 
=>- q + io 2 q + eq' 



q = p 



(3.3) 



q x +u 2 q 1 = -qo(t) 

Pi = ft 
q\ + u 2 q 2 = -3q 2 qi(t) 



P2 = <?2 



(3.4) 
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where q(t) = q cosut + — smut. The solutions for qi(t) and pi(t), are then given by 
(itsinut ( 2 tcosut £ 3 sin3u;£ (4COs3ujt 



2uo 2u 8u 2 8u 2 

Pi (t) = — (tuj cos uit + sin out) — — ( — too sin out + cos ut) — cos 3u;t H — — sin 3cut 

2ui 2uj 8uj 8lu 



(3.5) 



where the set of coefficients Ci, (2, (3, and £4 are given by 

t < \ 3 ( 3 , Q0P0 



4 V 4w 2 



4 V ^ 
l(pl 3qlp \ 



4 \w 3 4w / 



«*,«.) = j(^-flS) (3-6) 



As a consequence of the equations of motion, = 0. Additionally, a tedious but straight- 
forward exercise shows that {q(t),p(t)}p B = {q,p}pB- These observations imply that 
"time" t which appears in the solutions (3.5) plays the role of a parameter and that there 
is no physical significance to a particular instant in t. This justifies working with values of 
(q, p) at t = 0, without any loss of generality. Thus we will work at 0(e), with the solutions 

fl(r)\- n - $ 3 ^Po 
gi(0) - 9l -32^" 32^7 

Pl(0)=Pl = 32^ + W a7) 
and the corresponding ones at 0(e 2 ), which are given by 

= 5q o + 8q o p o + 17q ° p o 
92 256u 4 256a; 6 256cj 8 

= _56g^|_29g^o_ J7^_ 
P2 256a; 6 256u; 4 256a; 8 

We will proceed first with computations at the first order in perturbation theory. For 
notational simplicity, we introduce two parameters, (3 = an d « = f3ehu. Henceforth, 

// = y + —f 1 - Now using {q,p}pB = 1 and the solutions (3.7) we obtain for the Poisson 
bracket of the original variables 

Vpo} pB = (l-0eif o ) (3-9) 
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Thus the perturbation leads to a deformation of the Poisson algebra of the unperturbed 
observables (qo,po). This in particular implies that the Poisson bracket of arbitrary phase 
space functions f(qo,po) and g(qo,Po) is given by 

{f,g} pB = (l -/3e# ){/,s}p° } (3.10) 

where the superscript on the Poisson bracket implies zeroth order. With the help of the 
solutions the Hamiltonian at first order can be written as 

= H a +^HS (3.11) 
and at second order from (3.7) and (3.8) 

H(q ,Po) = H + + 2 -^fHl (3.12) 

We now turn to the quantum theory of this system. 



3.1 Symplectic structure, quantization and spectrum 

In order to obtain the spectrum of energies for the Hamiltonian (3.1 1) and (3.12), we need 
to promote the variables (qo,po) and the Hamiltonian H to hermitian operators on the 
Hilbert space H osc of the harmonic oscillator. The key to this quantization is the deforma- 
tion of the classical Poisson bracket algebra (3.9), which when promoted to a commutator 
via the correspondence, [, ] — > ih{, } yields 

[q ,Po] =th(l-f3eH ) (3.13) 

Observe that — 

1. The usual Heisenberg group HI in quantum mechanics is generated by the algebra 
[q,p\ = ih. Owing to the deformation of the classical Poisson algebra, q and p do 
not have this property. This renders the algebraic problem of operator ordering, such 
as Weyl ordering for example, unusually complex. 

2. The spirit of the present proposal is to use the Hilbert space of states defined by the 
unperturbed harmonic oscillator to construct the quantum theory of the perturbed 
problem at each order. Note, however, that covariant canonical quantization is the 
Heisenberg picture of quantum theory. It is well known that basis states of the 
Heisenberg picture evolve, and in this case, this evolution is imprinted as a defor- 
mation of the algebra. Hence, we need to to define a canonical pair generating a 
Heisenberg group at each order and to construct basis states \osc >G H QfiC , on which 
we can write the action of the Hamiltonian operator. 
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These observations motivate the transformation to a very simple canonical pair (q, p) that 
can be used to quantize the system at each order in perturbation theory. This transformation 
is what we referred to as an "adjustment" at the beginning of section 3. The new variables 
will be referred to throughout the remainder of this paper as a tilded-representation. The 
choice of this canonical pair is naturally motivated by the structure of the deformation 
of the classical Poisson algebra of the unperturbed oscillator variables. At 0(e) the new 
variables are 

q(qo,Po) : = qo^l + ^-H 

ft 
4 



p(q ,Po):=Po[l + ^H ) (3.14) 



This pair satisfies 



q,p\ =1 (3.15) 

PB 



The classical Hamiltonian (3.16) now becomes 



where 



H(q,p) : =H + ^H£ (3.16) 



Ho = ^ + ^f = H,-^Hl + 0^) (3.17) 



We define a quantum algebra by Weyl ordering the pair (3.14), 



Q ^ [q]w = Q= Ofo + Y [*>> H °\ +J 

P >-> [p]w = P= (po + Y [Po> H ] (3.18) 



To 0(e), therefore, we have a canonically conjugate pair which are functions of the oscil- 
lator data and which generates a Heisenberg group, 



q,p 



ih (3.19) 



A word in passing as to the meaning of Weyl ordering used in the above definition 
would be appropriate. Owing to the deformation in the commutator algebra eq. (3.9), 
we would expect a 0(e) modification of the standard Weyl ordering [29]. However, this 
does not contribute at 0(e) to the definitions (3.18), which incorporates the zeroth Weyl- 
ordered expression. Once these tilded-variables are defined, operator ordering to obtain 
0(e 2 ) expressions for operators corresponding to polynomial functions of (qo,po) does not 
pose a problem, since such functions can be re-expressed in the tilded representation by 
inverting the relations (3.14), using eq. (3.16) where appropriate. By construction the 
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operator of the tilded-representation generate a Heisenberg group, and hence usual Weyl 
ordering applies to functions of the form f(q m p r ). 
We now define a basis in terms of the operators 




UJ f z , 1 
uJ 




LO ( '* % 



2h 



[q p 

\ UJ 



Hn 1 



N = = -± - - (3.20) 

nuj 2 

where H = ^ + '^j-- Owing to eq.(3.19), [a u , at u ] = 1 and hence and at u satisfy the 
algebra of the usual creation-annihilation operators. As a consequence we can define states 

\n > that satisfy N\n >= h\n > with n E Z + . A straightforward calculation shows that 
the Hamiltonian (3.11) can be recast in terms of the Heisenberg pair (3.18), which upon 
yields the quantum Hamiltonian at 0(e), 

~ Be ~ o kHuj - 
H := Ho + ^-Hl + — r I (3.21) 
4 lb 

Given the definitions of eq. (3.20), the spectrum of H is given by (reinstating B and n) 

1\, 3en 2 /_o _ V 



H 



■n 



n H — 1 Huj H r- ( n 2 + n + 

2/ 8a; 2 V 2 



n) (3.22) 



The energy eigenvalues are in exact agreement with the results of Rayleigh-Schrodinger 
perturbation theory [25], implying thereby that our method reproduces the known result 
for the spectrum of the quartic oscillator in perturbation theory at this order. 

Computations to second order are straightforward, although algebraically more com- 
plex owing to the fact that the effect of the first order perturbation have to be taken into 
account at this order. First, the solutions (3.7) and (3.8) and eq. (3.10) imply the following 
deformation of the Poisson bracket of the unperturbed variables at second order: 

qo, Po ) =(l~ PeH + ^^H^) (3.23) 



pb V 64 

We observe already at this order that there is a emergent pattern in the way the perturbation 
appears as a deformation of the classical canonical algebra, the deformation being propor- 
tional to the unperturbed Hamiltonian raised to the power of the order of the perturbation. 
Given eq. (3.23) we can now define the tilded- pair at second order in analogy with (3.14) 
as 

Be TT 29p 2 e 2 Tj2 



*(1 + -4*0—3^0 
p:=Po(l + ^o-^# 2 ) (3-24) 
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We proceed with steps similar to those at first order. First, eq. (3.12) gives the full classical 
Hamiltonian at 0(e 2 ) , which when transformed to the tilded-representation becomes 



H = H + -H 



(3e ~ 2 17/3 2 e 2 ~ 



144 



m 



(3.25) 



Weyl ordering leads to the corresponding quantum Hamiltonian to 0(e 2 



H:=H + 



0e~ 



Hi + —I 



4 u 16 
which satisfies the eigenvalue equation, 



144 



-Hn H Hn 

576 



(3.26) 



H 



■n 



l\ t 3efr 2 / 2 _ 1 
n + - I nw + —— [n + n + - 

8u; 2 V 2 



e 2 h 3 
128cj 5 



(2n+ l) (^17n 2 + 17n + 21) 



(3.27) 



giving the quantized spectrum to second order that is once again in agreement with the 
standard quantum theory result. 

Having thus carried out computations to second order in perturbation theory, we are 
now in a position to suggest an algorithmic procedure for constructing a quantum theory 
to higher orders. 

To do so, observe that at first order the perturbation manifests as the deformation of 
the Poisson algebra of unperturbed variables (q , Po) given by eq. (3.9), which leads to the 
definition of the new variables given by eq. (3.14) which satisfy a canonical algebra at this 
order. Next, owing to eq. (3.10) and (3.23), if one now calculates the Poisson bracket of 
these new variables of first order, to include all terms of 0(e 2 ), one finds 



q,p 



1 



29/3 



2 £ 2 



PB 



64 



Hi 



(3.28) 



This deformation in the Poisson algebra of the new variables of first order leads to the 
definitions eq. (3.24) of the new variables at second order, where with a slight abuse of 
notation we rename variables that generate the Heisenberg group at each order to be (q, p). 
This suggests that for this model, we can iteratively define new variables that can be used 
to construct a quantum theory at each order in perturbation theory in the following way: 

If at any given order of perturbation n, we have a deformation of the commutative 
Poisson algebra, {Q, P} PB = (1 — «iA) of generic conjugate pairs Q = Q(qo,Po) and 
P = P(q ,Po), such that A(g , Po) is a homogeneous function of degree n, and ot\ ~ 0(e n ) 
is a constant — then, "new variables", 



q = Q 1 



p = P[l 



Oil 



n + 2 
n + 2 



-A 



-A 



(3.29) 
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satisfy a canonical algebra at that order. The corresponding ordered pair (q, p) generate a 
Heisenberg group. To see this note that — 



Pb 



Q,p) +^-(aq,p) (0) + ^-(q,ap' 

J pb n + 2{ J p B n + 2{ 

l-a 1 A > )+-^-(n + 2)A 
/ n + 2 



(3.30) 



where we have use the fact that at zeroth order Q = q and P = p and have exploited 

dA 
dqo 



Euler's theorem on homogeneous functions, which in particular means that g |A +po§^- 



nA. 

The conclusion is that at each order, with the help of the perturbative solutions a cor- 
responding quantum theory can be constructed and that with a choice of ordering (in this 
case Weyl ordering) we can pass at each level of perturbation theory from the classical 
evolving constants to their corresponding quantum counterparts. 

We now turn to a a discussion of the relevance of the results for quantum gravity. 



4 Conclusion and future outlook 

The results of the quartic oscillator model demonstrates through an explicit example that it 
is possible to extend covariant canonical quantization to include perturbation theory. Thus, 
I have demonstrated that where we do not a priori know the full classical solution, we can 
resort to a two-step procedure to find the corresponding quantum observables, namely 

1. Construct classical solutions via Hamilton's equations to perturbatively generate the 
solution space § at each order. 

2. Quantize the solution space thus generated using as the Hilbert space the time evolved 
basis states of the unperturbed problem. 

The result also addresses an issue raised in the past by Hajicek [24] regarding computabil- 
ity of "evolving constants", while at the same time providing a concrete example of how 
to implement of Bergmann's proposal for a covariant approximation scheme in quantum 
theory. An immediate application may be the quantum theory of minisuperspace models 
in general relativity which by virtue of their small number of degrees of freedom are quite 
tractable, yet physically relevant [30]. Recently such models have attracted a lot of at- 
tention as test beds of quantization techniques [31-34] . For instance, a problem of great 
physical interest is the study of the BKL behavior near cosmological singularities in quan- 
tum theory. Furthermore, dynamical system approaches where the Einstein equations can 
be reduced to a set of autonomous differential equations exist for such models [35,36]. 
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This opens up immediate possibilities of exploiting Hamiltonian techniques [26, 27] to 
study the quantum theory of such systems. If this bears fruit, then the next step would 
be to apply covariant canonical quantization to midisuperspaces, such as those of Gowdy 
models, whose classical data have recently been studied by several authors [37-41]. These 
questions are currently being investigated and will be reported in future publications. 
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